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1. Program Title : TOFAST - Fast Direct and Inverse Discrete Fourier Transform 

Routines. 

2. Abstract ; TOFAST is capable of calculating in place either the discrete 

- ^ vr 

Fourier transform (DFT) of real series jy.t , j = 0,1,...,N-1, where N = 2 
and M = 3,4,...,10, or the inverse discrete Fourier transform (DFT ) of the 
Fourier cosine and sine coefficients, a^, k = 0,1,...,n, and b^, k = 1,2, 
...,n-l, respectively, where n = N/2. Output of the DFT is the Fourier cosine 
and sine coefficients; output of the DFT ^ is the real series. Several 
modifications are possible which allow for input and/or output to be complex¬ 
valued . 


3. Hardware Requirement ; PDP-8 with minimal of 8K of core memory. (EAE is 
optional depending on version of floating-point package used (see Item 4)). 

4. Software Requirement ; One of the following extended memory floating-point 
packages must be in the memory field containing TOFAST, and must be able to 
be called by a JMS I 7: FPPES4, FPPES5, FPPES6, FPPEH4, FPPEH5, or FPPEH6. 
(The last three require EAE). 


5. Core Limits ; TOFAST uses 

Field 0 locations; 20-37, 64-77, 200-277, 400-1377; 

Field 1 locations: 200-(3*N+177), N is as in Item 2. 

Field 0 locations 25-37, 64-77 are used for variables and thus are available 

for that use by co-existing but independent programs defined by the user. 


6. Loading Procedure : 

(a) Load the extended memory floating-point package. 

(b) Load TOFAST binary tape. 


7. 


8 . 


Starting Addresses ; 

(a) For the DFT, SA is Field 0,0400. 

(b) For the DFT -1 , SA is Field 0,0401 


Use ; Prior to entry for either the DFT or the DFT \ the data array beginning 
at Field 1, 0200 must be properly filled: 

(a) For the DFT, the real series -£y^ is stored in j-sequence and without 

gaps, i.e, y0>y 1 »***,y N _] i * The data format can be either signed, 11-bit 
fixed-point, or standard floating-point. (Standard floating-point format 
refers to that described in the Floating-Point System Programming Manual, 
Digital 8-5-S). 


(b) For the DFT \ the Fourier cosine and sine coefficients are stored as the 


sequence a_ 


> a l* 


> a n * b l* b 2 **** ,b n-l* 1116 format is standard floating-point. 










With the data array filled as above, LOAD ADDress 000400 for the DFT, 
or 000401 for the DFT 1 . Set the switch register (SR) to N = size of data 
array; that is, for the DFT, N = length of series £y^ , and for the DFT - '*', 

N = total number of Fourier cosine and sine coefficients. In either case, 

M -1 

N must be of the form N = 2 , where M = 3,4,...,10. Press START. The DFT 

will run to completion, halting at 000551 with C(AC) = 0. The DFT, however, 
will halt at 000403 for format information. If the data format is floating¬ 
point, set SR to 4000; if fixed-point, set SR to number of bits to the right 
of the binary point. Press CONTinue. The DFT will run to completion and halt 
at 000444 with C(AC) =0. 

There is one error halt. It is at 000472 with C(AC) ?= 7777, and occurs 
if N is not of an acceptable form. Recovery is possible by resetting the 
SR to an acceptable N and pressing CONTinue. 


9. Operating Procedures : Assuming data array is filled according to Item 8, 

(a) for the DFT, 

(1) set SR to 000400; 

(2) Press LOAD ADDress; 

(3) set SR to N =no. of data in series; 

(4) press START; 

(5) at halt 000403,set SR to 
4000 if floating-point data, 

or * 

no. of bits right of binary point if fixed-point data; 

(6) press CONTinue: 

(b) for the DFT _1 , 

(1) set SR to 000401; 

(2) press LOAD ADDress; 

(3) set SR to N = total no. of cosine and sine coefficients; 

(4) press START. 


10 . 


Output : 
and sine 


b^jbgj.• 


At completion of the DFT, the data array contains the Fourier cosine 

» a l >•••> a n > 

.,b ,, and in standard floating-point format, where 

n-1 


coefficients, a^, 


and b^, respectively, in the sequence a^ 



- 2 - . 








and the series 


that was transformed can be represented by 


2~ 1 

yj - ¥ * I [% J k ) + \ si <^ J k )l + -§- cos (nj> 

k=0 

for j = 0 , X ,«•.)N“1• 

At completion of the DFT - *, the data array contains the series £y ^ f° r 
j = 0,1,...,N-l, in sequence and in standard floating-point format, and whose 
Fourier cosine and sine coefficients were contained in the input data array. 


11. Discussion : The method employed by TOFAST is a modification of the Cooley- 

Tukey fast Fourier transform algorithm. 

For the DFT, the series [y^ , 1= 0,1,...,N-l, and N = 2* 1 * M = 3,4, 

...,10, is first formed into the complex sequence 


V y aT + ly 2r-i 


for j = 0,1,...,n-l 


where n = N/2 = 2 , and where if the binary representation of j is 

J - then J 0 J l'" J n-2 J m-r LetCin8 xP den ° td che 

p-th intermediate vector, and defining Z(j,k) = j m _^j in _2*** 

I(J.k) - j n . 1 3 m .2-”Jk + l 1 Jk-l"-Jl'l0> and L <j’ k > " J -2 (">» d »). 


let 


and 


X? = X 


j for j = 0,1,...,n-l, 


x? = x ?- 1 


+ xr.: .. . w L (j* m -p ) 


i p-1 

Kj,p-D 


‘j ' Z(j,p-1) 

for j = ftl,...,n-l for every p = 1,2.. where W = e and i = JZ T. 

Then the complex Fourier coefficient at frequency* 0 ^ = 2llk/n is 

"ReCxJ] + Im[xJ] for k = 0 


Y <“k> ■ * ‘C / 1 + fr [5 C - <C ^ £or k - 1 - 2 . "- 1 




Re[X™] - Im[X™] for k - n 
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Finally 


XU 

where (X . ) is the complex congugate of X m . . 

n-K n-k 

\ = ~ReC Y (<t^)] for k = 0,1,...,n, 

and 

\ = - “Im[Y(c^)j for k = 1,2,...,n-1, 

^diere a k >b k are the Fourier cosine and sine coefficients, respectively, at 

*4 

frequency 2lIk/N. 

For the DFT , the Fourier cosine and sine coefficients, a and b, , res- 

k k’ 

pectively, are formed into the modified complex coefficients 


r 

for 

k = 0 

W - ib k 

for 

k = 1,2,...,n-l 

L 

L n 

for 

k = n 


where Y(aj^) = — Y^u^). Using the definitions and notations given for the 
DFT, let 


V k = 


■^[y^) + - Y (*« n )] for k = 0 

-k/2 

+ (Y(^ n _^>)*] - |r [Y(a, k > - (Y(w n k ))*] for k = 1,2, 


“ d Xj = Vy for j = 0,1,...,n-1 


. ,n-l 

# 


and 


X P = x p-i + x p_1 . w -L (J»m-p) 
S z<j,p-l) + X i(j,p-1) w 


for j “ 0>l»*>>»n-l for every p = 1,2, ... ,m. Finally, the series i 

given by 


y j = 


r Re[Xj /2 ] for j = 0,2,4,... ,N-2 
^ X (J-l)/2? for j = 1,3,5,...,N-1 
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12 . 


Some Possible User Modifications : It may be desired to locate TOFAST in 
field N and the data array in field M. To do this, the following changes to 
TOFAST should be made: 


REGISTER 

FROM 

TO 

CODING 

0441 

0017 

0N17 

FCDF N00 

0517 

0117 

0M17 

FCDF M00 

0546 

0017 

0N17 

FCDF N00 

0710 

6211 

62M1 

CDF M0 

0717 

6201 

62N1 

CDF N0 


In addition, several modifications to the extended memory floating-point 
package must be made in order to relocate it in field M. (See write-up for 
Overlay Modifications to the Floating-Point Packages.) 


There are four locations from which control can be transferred after the 
DFT and thfi DFT to a user-defined program; these four locations contain 
NOP's in the binary tapes supplied. After the DFT, jump from 0437 retaining 
the floating data field (FDF), or from 0443 not retaining the FDF; after the 
DFT 1 jump from 0544 retaining the FDF, or from 0550 not retaining the FDF. 

If the complex Fourier coefficients Y(a^) for the real series £y^ are 
desired as final output, make the following changes to TOFAST: 


REGISTER 

FROM 

TO 

CODING 

0430 

1255 

1353 

TAD LFNOP 

0432 

1255 

1353 

' TAD LFNOP 

0434 

1255 

1352 

TAD FGETP5 

0525 

1352 

1255 

TAD FDIVFN 

0527 

1352 

1255 

TAD FDIVFN 

0552 

3351 

5351 

FGETP5, 5351 

1234 

1026 

1030 

TAD NN 

1264 

5351 

5361 

FGET FN 

in Item 8, 

with N = length of 

I y jl • 

Output is the complex 


Fourier coefficients Y(ri^) = - iyf^, for k = 0,1, ...,n, stored as ^^, 

C^**..* 0t n *^l, <jr 2. ^i-I’ are transformed by the DFT -1 with the 

latter changes and using N = 2n, the output is the series £y ^ • 

If the series ^Y^ to be transformed is complex-valued make the following 
alterations to TOFAST: 
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REGISTER 

FROM 

TO 

CODING 

• 0412 

4616 

4617 

JMS I SPREAL 

0413 

4617 

4754 

JMS I CPXPML 

0422 

4620 

4754 

JMS I CPXPML 

0430 

1255 

5237 

JMP 437 

0505 

7240 

7000 

NOP 

0511 

7110 

7000 

NOP 

0522 

7770 

7774 

M4, -4 

0523 

6007 

7003 

M775, -775 

0525 

1352 

5334 

JMP534 

0543 

4620 

7000 

NOP 


13. 


Use 


fill data array with Yj - *, - ^ as ^ .S-l‘»#*»i.¥r — 

is as in Item 8 , with N = length of . Output is the Fourier coefficients 

Y(A^) = \ - i/^ k » stored as ot 0 » a 1 » • • • » ct n »^ 1 »^ 2 * * * * »/n-l* If Y (*k> are trans¬ 
formed by the DFT with the latter changes and using N = 2n, the output is 
the complex series stored as above. 

Logic Flow Charts : The two flow charts below are intended only to show the 
logic structures of the DFT and the DFT -1 of TOFAST; no attempt has been 

made to parallel in any way the physical structure of TOFAST and the flow 
charts. 


14 ‘ Program Listing : Listing follows flow charts. 

P-** References: The fast Fourier transform is discussed in: 

(1) Bingham, C., M.D. Godfrey, and J.W. Tukey, ‘'Modern Techniques of Power 
Spectrum Estimation." IEEE Trans. Audio Electroacousti cs. Vol. AU-15 
no. 2, (June 1967), pp. 56-66. 

(2) Bergland, G.D., "A Fast Fourier Transform Algorithm for Real-Valued 
Series." Commun. ACM , Vol. 11 , no.10, (Oct. 1968), pp. 703-710. 

(3) Cooley, J.W. and J.W. Tukey, "An Algorithm for the Machine Calculation 
of Fourier Series." Math. Comput .. Vol. 19, (1965), pp. 297-301. 

(4) Jennrich, R.I., B. Gilchrist, and A. Lynn, "REFAST: Fast Fourier Trans¬ 
form Routine." Program write-up of the Health Sciences Computing 
Facilities, UCLA, (Nov. 1967). 

(5) Rothman, J.E., "The Fast Fourier Transform and its Implementation." 

Decuscope , Vol. 7, no.3, (1968), pp. 3-10. > 

( 6 ) Singleton, R.C., "On Computing the Fast Fourier Transform." Commun. ACM. 

Vol. 10, no. 10, (Oct. 1967), pp. 646-654. - 

(7) Walton, P.L., "On the Fast Fourier Transform," in preparation. 
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Start of DFT; assumes data 
array (A) to be filled by 
real seriei ( Y i> such that 

A =Y. for \= 0 ),..., N-l, 

where N= 2 m =2 =2 m+ ^ / M is a 

n 

positive integer 



1.3 V 


SET: t = A. 

I 

V A ^ 

A k =t 


1.4 


LOOP 1 .0 

M3 ^ 


SET: A j=1 


FOR p* 1, 


m 


3.1 V 

SET: h = 

Ai 

^ i 

= 2(* j) 

f = 

7f/h 

t = 

cos f 

u = 

sin f 

c = 

1 

s = 

0 


TOFAST LOGIC FLOW CHART 
PAGE 1 OF 6 


j is j bit-reversed as an 
m-bit integer, i.e. if j=j m 


' m - 2 "V 


then 
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i 


! 


i 
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T OF AST LOGIC 
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FOR p: 1/m 

5 


5.1 


SET: h = 4 j 

A i = 2 < 4 i) 
f= ir/h 

t = cos f 
u = sin f 
c = 1 

s=0 


5.20 


5.210 


FOR I: = 0/h-1 


FOR j: = l,n-l;step A j 
1 ? 


5.212 


LOOP 5.210 


5.22 y 


sTHT 

C 

S 

= c 

= r.t-s.u 

= s. Hr. u 




© 


TOFAST LOGIC FLOW CHART 
PAGE 5 OF 6 


SET: k = j + h 

d = c.A k -s.A n+k 

e = C * A n + k + S - A k 

A, = A. - d 
k 1 

A n + k = A n +i - e 

A. = A. + d 

1 1 


Note sign changes in 

calculating d, e for 
.DFT and DFT this 

is major difference^ 
between transformaron 
subroutines 



A = A . + e 
n+| n+| 



± 
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/PROGRAM TOFAST - FAST FORWARD OR INVERSE DISCRETE 
/FOURIER TRANSFORM FOR REAL DATA 
/PAGE 0 OF TOFAST 
*20 

/TOFAST CONSTANTS 

/SA OF DATA ARRAY 


0020 0200 

DATAST, 

200 

0021 0002 

PI/ 

2 

0022 3110 


3110 

0023 3755 


3755 

0024 0262 

INDEXL, 

INDEX 


/TOFAST VARIABLES 

0025 0000 

M, 

0 

0026 0000 

N, 

0 

002 7 0000 

NMIN1, 

0 

0030 0000 

NN, 

0 

0031 0000 

J, 

0 

0032 0000 

K, 

0 

0033 0000 

AJL, 

0 

0034 0000 

ANPJL, 

0 

0035 0000 

AKL, 

0 

0036 0000 

ANPKL, 

0 

0037 0000 

H, 

0 

*64 

0064 0000 

c. 

0 

0065 0000 


0 

0066 0000 


0 

0067 0000 

s. 

0 

0070 0000 


0 

0071 0000 


0 

0072 0000 

T/ 

0 

0073 0000 


0 

0074 0000 


0 

0075 0000 

u. 

0 

0076 0000 


0 

0077 0000 


0 


/PAGE 1 OF TOFAST 
*200 


/SUBROUTINE TO DO PARTIAL 
/VECTOR 

0200 0000 

LOOP3, 

0 

0201 1031 


TAD J 

0202 1037 


TAD H 

0203 3032 


DCA K 

0204 4242 


JMS SETLOC 

0205 4407 


FENTER 

0206 5436 


FGET 1 ANPKL 

0207 3067 


FMPYS 

0210 6272 


FPUT D1 

0211 5435 


FGET 1 AKL 

0212 3064 


FMPYC 


/K=JfH 

/GET ADDRESSES 
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,0213 1272 

SWIT1, 

FADD D1 

/OR FSUB D1 'FOR INVERSE 

0214 6272 


FPUT D1 

/D1=C*A(K)+S*A(N+K) 

0215 5435 


FGET 1 AKL 

/OR D1=C*A(K)-S*A(N+K) 

0216 3067 


FMPYS 


0217 6275 


FPUT El 


0220 5436 


FGET 1 ANPKL 


0221 3064 


FMPYC 


0222 2275 

SWIT2, 

FSUB El 

/OR FADD El FOR INVERSE 

0223 6275 


FPUT El 

/El=C*A(NfK)-S*A(K) 

0224 5433 


FGET 1 AJL 

/OR E1=C*A(N+ K)+S*A(K) 

0225 2272 


FSUB D1 


0226 6435 


FPUT 1 AKL 

/A(K)=A(J)-D1 

0227 5434 


FGET 1 ANPJL 


0230 2275 


FSUB El 


0231 6436 


FPUT 1 ANPKL 

/A(NfK)=A(NfJ)-El 

0232 5433 


FGET 1 AJL 


0233 1272 


FADD D1 


0234 6433 


FPUT 1 AJL 

/A(J)=A(JHD1 

0235 5434 


FGET 1 ANPJL 


0236 1275 


FADD El 


0237 6434 


FPUT 1 ANPJL 

/A(Nf J)=A(N+J)+E1 

0240 0000 


FEXT 


0241 5600 


JMP 1 LOOP3 

/EXIT 


/SUBROUTINE TO FIND FLOATING ADDRESSES OF A(J),A(NfJ), 


/A(K), AND A(N+K) 


0242 0000 

SETLOC, 

0 


0243 1031 


TAD J 


0244 4262 


JMSINDEX 


0245 3033 


DCA AJL 

/ADDRESS OF A(J) 

0246 1026 


TAD N 


0247 1031 


TAD J 


0250 4262 


JMS INDEX 


0251 3034 


DCA ANPJL 

/ADDRESS OF A(N+J) 

0252 1032 


TAD K 


0253 4262 


JMSINDEX 


0254 3035 


DCA AKL 

/ADDRESS OF A(K) 

0255 1026 


TAD N 


0256 1032 


TAD K 


0257 4262 


JMS INDEX 


0260 3036 


DCA ANPKL 

/ADDRESS OF A(N+K) 

0261 5642 


JMP 1 SETLOC 

/EXIT 


/SUBROUTINE TO FIND FLOATING ADDRESS OF AtC(AC)3 

0262 0000 

INDEX, 

0 


0263 3271 


DCA INXSTR 


0264 1271 


TAD INXSTR 


0265 7104 


CLL RAL 


0266 1271 


TAD INXSTR 


0267 1020 


TAD DATAST 


0270 5662 


JMP 1 INDEX 
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0271 0000 

INXSTR, 

0 


0272 0000 

Dl, 

0 


0273 0000 


0 


0274 0000 


0 

\ 

0275 0000 

El, 

0 


0276 0000 


0 


0277 0000 


0 



/PAGE 2 OF TOFAST 




*400 



/CONTROL ROUTINE; ENTER AT 0400 FOR DFT; 


/ENTER AT 0401 FOR INVERSE DFT 


0400 5202 


JMP FWDFFT 


0401 5324 


JMP REVFFT 



/DFT; ACCEPTS 2TMP1 FIXED- OR FLOATING-POINT DATA, 


/WHERE MP1 

=3,4,5,...,10 


0402 4257 

FWDFFT, 

JMS PRELIM 

/SET UP CONSTANTS, FIELD 

0403 7402 


HLT 

/HALT FOR DATA FORMAT 

0404 7604 


LAS 

/BIT 0 OF SR SET IF 

0405 7510 


SPA 

/FLOATING;OTHERWISE SET 

0406 5221 


JMP DATFLT 

/TO # BITS RIGHT OF . 

0407 7041 


CIA 

/FIXED-POINT DATA 

0410 1215 


TAD Cl 3 

/#BITS LEFT OF . 

0411 3757 


DCA 1 BINPNL 

0412 4616 


JMS 1 FIXPML 

/PERMUTE TO COMPLEX 

0413 4617 


JMS 1 SPREAL 

/VECTOR,THEN FLOAT 

0414 5223 


JMP DODFT 


0415 0013 

Cl 3, 

13 


0416 0600 

FIXPML, 

FIXPMT 


0417 1121 

SPREAL, 

SPREAD 


0420 0606 

FLTPML, 

FLTPMT 


0421 7300 

DATFLT, 

CLA CLL 

/FLOATING-POINT DATA 

0422 4620 


JMS 1 FLTPML 

/PERMUTE TO COMPLEX 

0423 1245 

DODFT, 

TAD FADDD1 


0424 3647 


DCA 1 SWIT1L 


0425 1246 


TAD FSUBE1 


0426 3650 


DCA 1 SWIT2L 

/SET LOGIC FOR DFT 

0427 4654 


JMS 1 MAINL 

/DO DFT 

0430 1255 


TAD FDIVFN 

J 

0431 3651 


DCA 1 SWIT3L 


0432 1255 


TAD FDIVFN 


0433 3652 


DCA 1 SWIT4L 


0434 1255 


TAD FDIVFN 


0435 3653 


DCA 1 SWIT5L 

/SET LOGIC FOR UNSCRAM. 

0436 4656 


JMS 1 UNSCRL 

/UNSCRAMBLE 

0437 7000 


NOP 

/JUMP FROM HERE TO USER- 

0440 4407 


FENTER 

/DEFINED ROUTINES 

0441 0017 


FCDF0 

/SET FLOATING DATA 

0442 0000 


FEXT 

/FIELD=INST. FIELD 

0443 7000 


NOP 

/JUMP FROM HERE, TOO 

0444 7402 


HLT 

/HALT 
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0445 1272 

FADDD1, 

1272 

/FADD D1 ON PAGE 1 

0446 2275 

FSUBE1, 

2275 

/FSUB El ON PAGE 1 

0447 0213 

SWIT1L, 

SWIT1 


0450 0222 

SWIT2L, 

SWIT2 


0451 1244 

SWIT3L, 

SWIT3 


0452 1250 

SWIT4L, 

SWIT4 


0453 1265 

SWIT5L, 

SWIT5 


0454 1000 

MAINL, 

MAIN 


0455 4361 

FDIVFN, 

4361 

/FDIV FN ON PAGE 5 

0456 1200 

UNSCRL, 

UNSCRM 



/SUBROUTINE TO SET UP CONSTANTS FOR TRANSFORM, 


/AND SET FLOATING DATA FIELD 


0457 0000 

PRELIM, 

0 


0460 7604 

LAS 


/ENTERING SR CONTAINS 

0461 3030 


DCA NN 

/NN=NO. of real data 

0462 1030 


TAD NN 


0463 1322 


TAD M10 


0464 7510 


SPA 


0465 5271 


JMP ERRNN 


0466 1323 


TAD Ml 771 


0467 7710 


SPA CLA 

/IS 7<NN<1025 

0470 5274 


JMP .+4 


0471 7240 

ERRNN, 

STA 

/NO-HALT WITH 7777 IN 

0472 7402 


HLT 

/AC; RESET SR TO LEGAL 

0473 5260 


JMP PRELIM+1 

/NN, PRESS CONTINUE 

0474 3760 


DCA 1 MP1L 

/YES 

0475 1030 


TAD NN 


0476 7110 


CLL RAR 


0477 7430 


SZL 


0500 5303 


JMP .+3 


0501 2760 


ISZ 1 MP1L 


0502 5276 


JMP .-4 


0503 7640 


SZA CLA 

/IS NN OF FORM 2fMPl 

0504 5271 


JMP ERRNN 

/NO-ERROR HALT 

0505 7240 


STA 

/YES 

0506 1760 


TAD 1 MP1L 


0507 3025 


DCA M 

/M=MP1-1 

0510 1030 


TAD NN 


0511 7110 


CLL RAR 


0512 3026 


DCA N 

/N=Nls|/2 

0513 7240 


STA 


0514 1026 


TAD N 


0515 3027 


DCA NMIN1 

/N-l 

0516 4407 


FENTER 


0517 0117 


FCDF 100 

/SET FOR FLOATING DATA 

0520 0000 


FEXT 

/FIELD=1 

0521 5657 


JMP 1 PRELIM 

/EXIT 
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0522 7770 
0523 6007 


0524 4257 
0525 1352 
0526 3651 
0527 1352 
0530 3652 
0531 1353 
0532 3653 
0533 4656 
0534 4754 
0535 1355 
0536 3647 
0537 1356 
0540 3650 
0541 4654 
0542 4754 
0543 4620 
0544 7000 
0545 4407 
0546 0017 
0547 0000 
0550 7000 
0551 7402 
0552 3351 
0553 0010 
0554 0625 
0555 2272 
0556 1275 
0557 1176 
0560 0765 


0600 0000 
0601 1205 
0602 3253 
0603 4214 
0604 5600 
0605 4301 

0606 0000 
0607 1213 
0610 3253 
0611 4214 
0612 5606 


M10, -10 

Ml 771, -1771 

/INVERSE DFT; ACCEPTS 2tMPl FLOATING-POINT FOURIER 
/COEFFICIENTS, WHERE MP1=3,4,5,..10, AND COEFF.S 


REVFFT, JMS PRELIM 

/SET UP CONSTANTS, FIELD 

TAD FMPYP5 

DCA 1 SWIT3L 
TAD FMPYP5 

DCA 1 SWIT4L 
TAD LFNOP 

DCA 1 SWIT5L 

/SET LOGIC FOR SCRAMBLE 

JMS 1 UNSCRL 

/SCRAMBLE 

JMS 1 CPXPML 

/PERMUTE COMPLEX VECTOR 

TAD FSUBD1 

DCA 1 SWIT1L 
TAD FADDE1 

DCA 1 SWIT2L 

/SET LOGIC FOR INVERSE 

JMS 1 MAINL 

/DO INVERSE DFT 

JMS 1 CPXPML 
JMS 1 FLTPML 

/PUT DATA IN- SEQUENCE 

NOP 

/JUMP FROM HERE TO USER- 

FENTER 

/DEFINED ROUTINES 

FCDF0 

/SET FLOATING DATA 

FEXT 

/field=inst. field 

NOP 

/OR FROM HERE 

HLT 

/HALT 

FMPYP5, 3351 

/FMPY ZP5 ON PAGE 5 

LFNOP, FNOP 

CPXPML, CPXPMT 

FSUBD1, 2272 

/FSUB D1 ON PAGE 1 

FADDE1, 1275 

/FADD El ON PAGE 1 

BINPNL, BINPNT 

MP1L, MP1 

/PAGE 3 OF TOFAST 
*600 

/SUBROUTINE TO BIT-REVERSE 

PERMUTE NN FIXED-POINT DATA 

FIXPMT, 0 

TAD JMSFXP 

DCA PMTMOV 

/SET LOGIC 

JMS PMTDAT 

/PERMUTE 

JMP 1 FIXPMT 

/EXIT 

JMSFXP, JMS FIXMOV 

/SUBROUTINE TO BIT-REVERSE 

PERMUTE NN FLOATING DATA 

FLTPMT, 0 

TAD JMSFLT 

DCA PMTMOV 

/SET LOGIC 

JMS PMTDAT 

/PERMUTE 

JMP 1 FLTPMT 

/EXIT 
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0613 

4321 

JMSFLT, 

JMS FLTMOV 




/SUBROUTINE TO SET UP FOR PERMUTION OF NN REAL DATA 

0614 

0000 

PMTDAT, 

0 


0615 

1224 


TAD LNOP 


0616 

3275 


DCA FLPCTL 

/SET LOGIC 

0617 

1027 


TAD NMIN1 


0620 

7104 


CLL RAL 


0621 

7041 


CIA 

/NN-2 DATA 

0622 

4240 


JMS PERMUT 

/PERMUTE 

0623 

5614 


JMP 1 PMTDAT 

/EXIT 

0624 

7000 

LNOP, 

NOP 




/SUBROUTINE TO BIT-REVERSE PERMUTE N FLOATING-POINT 



/COMPLEX DATA 


0625 

0000 

CPXPMT, 

0 


0626 

1237 


TAD JMSCPX 


0627 

3253 


DCA PMTMOV 


0630 

1267 


TAD CLLRAR 


0631 

3275 


DCA FLPCTL 

/SET LOGIC 

0632 

7240 


STA 


0633 

1027 


TAD NMIN1 


0634 

7041 


CIA 

/N-2 DATA 

0635 

4240 


JMS PERMUT 

/PERMUTE 

0636 

5625 


JMP 1 CPXPMT 

/EXIT 

0637 

4341 

JMSCPX, 

JMS CPXMOV 




/SUBROUTINE TO BIT-REVERSE PERMUTE 

-C(AC) DATA AS 



/SET UP BY LOGIC ALTERATIONS 


0640 

0000 

PERMUT, 

0 


0641 

3360 


DCA CTR1 

/-NO. OF DATA 

0642 

7201 


CLA IAC 


0643 

3031 


DCA J 

/START WITH >1 

0644 

1031 

PMTNXT, 

TAD J 


0645 

4260 


JMS FLIP 


0646 

3032 


DCA K 

/K=BIT-REVERSED j 

0647 

1032 


TAD K 


0650 

7041 


CIA 


0651 

1031 


TAD J 


0652 

7710 


SPA CLA 

/IS K>J 

0653 

0000 

PMTMOV, 

0 

/YES-SWITCH DATA J,K 

0654 

2031 


ISZ J 

/NEXT J 

0655 

2360 


ISZ CTR1 

/HAVE ALL DATA BEEN DONE 

0656 

5244 


JMP PMTNXT 

/NO-NEXT 

0657 

5640 


JMP 1 PERMUT 

/YES-EX IT 



/SUBROUTINE TO BIT REVERSE C(AC) AS MP1- OR M-BIT WORD 

0660 

0000 

FLIP, 

0 


0661 

3277 


DCA FLPJ 


0662 

1365 


TAD MP1 


0663 

7041 


CIA 


0664 

3361 


DCA CTR2 

/MP1 BITS TO BE REV. 
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0665 3300 
0666 1277 
0667 7110 
0670 3277 
0671 1300 
0672 7004 
0673 2361 
0674 5265 
0675 7000 
0676 5660 
0677 0000 
0700 0000 


FLPLOP, 

DCA FLPK 

TAD FLPJ 


CLLRAR, 

CLL RAR 

DCA FLPJ 

TAD FLPK 

/READ J FROM R TO L 


RAL 

/BUILD K FROM L TO R 


ISZ CTR2 

/HAVE ALL BITS BEEN REV. 

FLPCTL, 

JMP FLPLOP 

/NO-NEXT 

NOP 

/YES (CLL RAR FOR M-BIT) 


JMP 1 FLIP 

/EXIT 

FLPJ, 

0 


FLPK, 


0 


/SUBROUTINE TO SWITCH FIXED-POINT A(J) AND A(K) 


0701 0000 FIXMOV, 0 


1031 

TAD J 


1020 

TAD DATAST 


3033 

DCA AJL 

/ADDRESS OF A(J) 

1032 

TAD K 


1020 

TAD DATAST 


3035 

DCA AKL 

/ADDRESS OF A(K) 

6211 

CDF 10 


1435 

TAD 1 AKL 


3362 

DCA T3 

/T3=A(K) 

1433 

TAD 1 AJL 


3435 

DCA 1 AKL 

/A(K)=A(J) 

1362 

TAD T3 


3433 

DCA 1 AJL 

/A(J)=T3 

6201 

CDF 0 


5701 

JMP 1 FIXMOV 

/EXIT 


0721 0000 
0722 1031 
0723 4424 
0724 3033 
0725 1032 
0726 4424 
0727 3035 
0730 4407 
0731 5435 
0732 6362 
0733 5433 
0734 6435 
0735 5362 
0736 6433 
0737 0000 
0740 5721 


FLTMOV, 0 


TAD J 

JMS 1 INDEXL 

DCA AJL 

TAD K 

/ADDRESS OF A(J) 

JMS 1 INDEXL 

DCA AKL 

FENTER 

/ADDRESS OF A(K) 

FGET 1 AKL 

FPUT T3 

FGET 1 AJL 

/T3=A(K) 

FPUT 1 AKL 

FGET T3 

/A(K)=A(J) 

FPUT 1 AJL 

FEXT 

/A(J)=T3 

JMP 1 FLTMOV 

/EXIT 


/SUBROUTINE TO SWITCH COMPLEX DATA S(J) AND S(K) 
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0741 

0000 

CPXMOV, 

0 


0742 

4321 


JMS 

FLTMOV 

0743 

1031 


TAD 

J 

0744 

3357 


DCA 

SAVEJ 

0745 

1026 


TAD 

N 

0746 

1031 


TAD 

J 

0747 

3031 


DCA 

J 

0750 

1026 


TAD 

N 

0751 

1032 


TAD 

K 

0752 

3032 


DCA 

K 

0753 

4321 


JMS 

FLTMOV 

0754 

1357 


TAD 

SAVEJ 

0755 

3031 


DCA 

J 

0756 

5741 


JMP 

1 CPXMOV 

0757 

0000 

SAVEJ, 

0 


0760 

0000 

CTR1, 

0 


0761 

0000 

CTR2, 

0 


0762 

0000 

T3, 

0 


0763 

0000 


0 


0764 

0000 


0 


0765 

0000 

MP1, 

0 




/PAGE 4 OF 

TOFAST 


/SWITCH REAL PARTS 
/SAVE J 


/IM OF S(J)=A(NH-J) 


/IM OF S(K)=A(NfK) 
/SWITCH IMAGIN. PARTS 

/RESTORE J 
/EXIT 


*1000 

/SUBROUTINE TO TRANSFORM N FLOATING-POINT COMPLEX DATA 


1000 0000 

MAIN, 

0 


1001 7201 


CLA IAC 


1002 3367 


DCA DJ 

/START WITH DJ=1 

1003 1025 


TAD M 


1004 7041 


CIA 


1005 3370 


DCA MINP 

/M INTERMEDIATE TRANS. 

1006 4236 

PLOOP, 

JMS LOOP1 

/SET UP FOR 1 NT. TRANS. 

1007 3371 


DCA L 

/L=0 

1010 1037 


TAD H 


1011 7041 


CIA 


1012 3372 


DCA MINL 

/H PART./INT. 

1013 1371 

LLOOP, 

TAD L 


1014 3031 


DCA J 

/ 

1015 4635 

JLOOP, 

JMS 1 LOOP3L 

/DO INT. PARTIAL 

1016 1031 


TAD J 


1017 1367 


TAD DJ 


1020 3031 


DCA J 

/J=>DJ 

1021 1031 


TAD J 


1022 7041 


CIA 


1023 1027 


TAD NMIN1 


1024 7700 


SMA CLA 

/IS J>N-1 

1025 5215 


JMP JLOOP 

/NO-NEXT PASS OF PART. 

1026 2372 


ISZ MINL 

/YES-ARE ALL PART.S DONE 

1027 5233 


JMP MOREL 

/NO-NEXT PARTIAL 
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1030 2370 

1031 5206 

1032 5600 

1033 4273 

1034 5213 

1035 0200 

1036 0000 

1037 1367 

1040 3037 

1041 1367 

1042 7104 

1043 3367 

1044 1037 

1045 3045 
i 1046 4407 

1047 0011 

1050 6373 

1051 5021 

1052 4373 

1053 6373 

1054 0004 

1055 6072 

1056 5373 

1057 0003 

1060 6075 

1061 0000 

1062 7201 

1063 3064 

1064 7132 

1065 3065 

1066 3066 
C 1067 3067 

1070 3070 

1071 3071 

1072 5636 


1073 0000 

1074 4407 

1075 5064 

1076 6373 

1077 5067 

1100 3075 

1101 6064 

1102 5373 

1103 3072 


ISZ MINP 
JMP PLOOP 
JMP I MAIN 
MOREL, JMS SETUPL 
JMP LLOOP 
LOOP3L, LOOP3 
/SUBROUTINE TO SET UP 
LOOP1, 0 


/YES-ARE ALL INT.S DONE 
/NO-NEXT INTERMEDIATE 
/YES-EX IT 

/SET UP FOR NEXT PARTIAL 


FOR INTERMEDIATE TRANSFORM 


TAD DJ 

DCA H /H=DJ 

TAD DJ 
CLL RAL 

DCA DJ /DJ=2*DJ 

TAD H 

DCA 45 

F ENTER 

FLOAT 

FPUT R4 

FGET PI 

FDIV R4 

FPUT R4 /R4=PI/H 

FCOS 

FPUT T /T=CO S(R4) 

FGET R4 

FSIN 

FPUT U /U=SIN(R4) 

FEXT 

CLA IAC 

DCA C 

STL RTR 

DCA C+l 

DCA C+2 /c=T 

DCA S 
DCA S+1 
DCA S+2 


JMP I LOOP1 /EXIT 

/SUBROUTINE TO SET UP FOR INTERMEDIATE PARTIAL TRANS¬ 
FORM EXCEPT FOR FIRST ONE 
SETUPL, 0 


FENTER 
FGET C 

FPUT R4 /R4=C 

FGET S 

FMPY U 

FPUT C 

FGET R4 

FMPY T 
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1105 6064 


FPUT C 

/C=R4*T-S*U 

1106 5067 


FGET S 


1107 3072 


FMPYT 


1110 6067 


FPUT S 


1111 5373 


FGET P4 


1112 3075 


FMPY U 


1113 1067 


FADD S 


1114 6067 


FPUT S 

/S=S*T+R4*U 

1115 0000 


FEXT 


1116 2371 


ISZ L 

/L=L+1 

1117 7000 


NOP 


1120 5673 


JMP 1 SETUPL 

/EXIT 


/SUBROUTINE TO FLOAT NN FIXED- 

■POINT DATA 

1121 0000 

SPREAD, 

0 


1122 1030 


TAD NN 


1123 7041 


CIA 


1124 3377 


DCA MINJ4 

/NN DATA 

1125 7240 


STA 


1126 1030 


TAD NN 


1127 3031 


DCA J 

/START AT END 

1130 1031 


TAD J 


1131 4424 


JMS 1 INDEXL 


1132 3365 


DCA SPDFLT 

/ADD. OF LAST FLOAT 

1133 1031 


TAD J 


1134 1020 


TAD DATAST 


1135 3366 


DCA SPDFIX 

/ADD. OF LAST FIXED 

1136 6211 

SPDNXT, 

CDF 10 


1137 1766 


TAD 1 SPDFIX 


1140 6201 


CDF 0 


1141 3045 


DCA 45 


1142 3046 


DCA 46 


1143 3047 


DCA 47 


1144 1376 


TAD BINPNT 


1145 3044 


DCA 44 


1146 4764 


JMS1DNORML 

/FLOAT DATA 

1147 4407 


FENTER 


1150 6765 


FPUT 1 SPDFLT 

/PUT IT AWAY 

1151 0000 


FEXT 


1152 7240 


STA 


1153 1366 


TAD SPDFIX 


1154 3366 


DCA SPDFIX 

/RETREAT POINTERS 

1155 7126 


STL RTL 


1156 7040 


CMA 


1157 1365 


TAD SPDFLT 


1160 3365 


DCA SPDFLT 


1161 2377 


ISZ MINJ4 

/ARE ALL FLOATING 

1162 5336 


JMP SPDNXT 

/NO-NEXT 

1163 5721 


JMP 1 SPREAD 

/YES-SWIM BACK 





1164 6600 

1165 0000 

1166 0000 
1167 0000 

1170 0000 

1171 0000 

1172 0000 

1173 0000 

1174 0000 

1175 0000 

1176 0000 

1177 0000 


DNORML, 
SPDFLT, 
SPDFIX, 
DJ, 

MINP, 

L/ 

MINL, 

R4/ 


DNORM 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

/PAGE 5 OF TOFAST 

' 1 200 


BINPNT, 

MINJ4, 


1200 0000 

/SUBROUTINE TO UN- OR RE- 
UNSCRM, 0 

1201 4230 


JMS SETUPF 

1202 1026 


TAD N 

1203 7110 


CLL RAR 

1204 7041 


CIA 

1205 7001 


IAC 

1206 3037 


DCA H 

1207 7201 


CLA IAC 

1210 3031 


DCA J 

1211 4271 

UNLOOP, 

JMS LOOP4 

1212 2031 

ISZ J 

1213 2037 


ISZ H 

1214 5211 


JMPUNLOOP 

1215 2361 


ISZ FN 

1216 4755 


JMS 1 SETLOL 

1217 4407 


FENTER 

1220 5433 


FGET 1 AJL 

1221 3361 


FMPY FN 

1222 6433 


FPUT 1 AJL 

1223 5434 


FGET 1 ANPJL 

1224 3361 


FMPY FN 

1225 6434 


FPUT 1 ANPJL 

1226 0000 


FEXT 

1227 5600 


JMP 1 UNSCRM 


•SCRAMBLE COEFFICIENTS 

/DO A(0),A(N);SET UP REST 


/N/2-1 ITERATIONS 

/START WITH J= 1 
/UN- OR RE-SCRAMBLE 
/NEXT J 

/IS SCRAMBLING DONE 
/NO-NEXT 

/YES-DO A(N/2),A(N+N/2) 


/SET UP FOR REST 

1230 0000 SETUPF, 0 

T231 1,026 TAD N 

1232 4424 JMS I INDEXL 

1233 3035 DCA AKL 

1234 1026 TAD N 

1235 3045 DCA 45 

1236 4407 FENTER 


/EXIT 

SCRAMBLE A(0),A(N), AND 


/ADDRESS OF A(N) 
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1237 0011 


FLOAT 

1240 6361 


FPUT FN 

1241 5420 


FGET 1 DATAST 

1242 6072 


FPUT T 

1243 1435 


FADD 1 AKL 

1244 4361 

SWIT3, 

FDIV FN 

1245 6420 


FPUT 1 DATAST 

1246 5072 


FGET T 

1247 2435 


FSUB 1 AKL 

1250 4361 

SWIT4, 

FDIV FN 

1251 6435 


FPUT 1 AKL 

1252 5021 


FGET PI 

1253 4361 


FDIV FN 

1254 6372 


FPUT F5 

1255 0004 


FCOS 

1256 6072 


FPUT T 

1257 6064 


FPUT C 

1260 5372 


FGET F5 

1261 0003 


FSIN 

1262 6075 


FPUT U 

1263 6067 


FPUT S 

1264 5351 


FGET ZP5 

1265 4361 

SWIT5, 

FDIV FN 

1266 6361 


FPUT FN 

1267 0000 


FEXT 

1270 5630 


JMP 1 SETUPF 


/SUBROUTINE TO DO ACTU/ 

1271 0000 

LOOP4, 

0 

1272 1031 


TAD J 

1273 7041 


CIA 

1274 1026 


TAD N 

1275 3032 


DCA K 

1276 4755 


JMS 1 SETLOL 

1277 4407 


FENTER 

1300 5434 


FGET 1 ANPJL 

1301 1436 


FADD 1 ANPKL 

1302 6375 


FPUT R5 

1303 3064 


FMPY C 

1304 6364 


FPUT D5 

1305 5375 


FGET R5 

1306 3067 


FMPY S 

1307 6375 


FPUT R5 

1310 5433 


FGET 1 AJL 

1311 2435 


FSUB 1 AKL 

1312 6356 


FPUT O 

1313 3064 


FMPY C 

1314 1375 


FADD R5 

1315 6367 


FPUT E5 


/T=A(0) 

/OR FMPY ZP5 FOR INVERSE 
/A(0)<A(0)+A(N)3/N / 

/OR A(0)=C A(0)+A(N) 3/2 

/OR FMPY ZP5 FOR INVERSE 
/A(N)= C T-A(N) 3/N, 

/OR A(N)= CT T-A(N) 3/2 

/F5=PI/N 

/T=CO S(F5) 

/ C=T 


/U=SIN(F5) 

/S=U 

/OR FNOP FOR INVERSE 
/FN=.5/N, OR FN=.5 

/EXIT 

I- OR RE-SCRAMBLING 


/K=N-J 

/GET ADDRESSES 


/R5=A(N+J)+A(N+K) 

/D5=C*R5 


/R5=S*R5 


/Q=A(J)-A(K) 


/E5=C*Q+R5 


i 
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1316 5356 

1317 3,067 

1320 2364 

1321 6356 

1322 5433 

1323 1435 

1324 6435 

1325 2356 

1326 3361 

1327 6433 

1 330 5435 

1331 1356 

1332 3361 

1333 6435 

1334 5434 

1335 2436 

1336 6434 

1337 1367 

1340 3361 

1341 6436 

1342 5367 

1343 2434 

1344 3361 

1345 6434 


1346 0000 

1347 4754 
135^ 5671 


1351 0000 ZP 5, 

1352 2000 

1353 0000 


1354 

1073 

SETULL, 

1355 

0242 

SETLOL, 

1356 

0000 

Q, 

1357 

0000 


1360 

0000 


1361 

0000 

FN, 

1362 

0000 


1363 

0000 


1364 

0000 

D5, 

1365 

0000 


1366 

0000 


1367 

0000 

E5, 

137J0 

0000 


1371 

0000 


1372 

0000 

F 5, 

1373 

0000 


1374 

0000 



FGET Q 
FMPYS 
FSUB D5 
FPUTQ 
FGET I AJL 
FADD I AKL 
FPUT I AKL 
FSUB Q 
FMPY FN 
FPUT I AJL 
FGET I AKL 
FADD Q 
FMPY FN 
FPUT I AKL 
FGET I ANPJL 
FSUB I ANPKL 
FPUT I ANPJL 
FADD E5 
FMPY FN 
FPUT I ANPKL 
FGET E5 
FSUB I ANPJL 
FMPY FN 
FPUT I ANPJL 
FEXT 

JMS I SETULL 
JMP I LOOP4 
0 

2000 

0 

SETUPL 

SETLOC 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


/ Q=S*Q-D5 
/A(K)=A(J)+A(K) 

/A(J)=Ca(K)-Q3*FN 

/A(K)=D A(K)+Q 3 *FN 
/A(N+J)=A(NfJ)-A(N+K) 
/A(N+K)=C A(NfJ)+E5 3 *fn 

/A(N+J)=C E5-A(Nf J) 3*FN 

/SET UP FOR NEXT GROUP 
/EXIT 
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1375 0000 R5, 0 

1376 0000 0 

1377 0000 0 

/DEFINITIONS FOR TOFAST 
FENTER=JMS I 7 
D NORM-66/0 

$ 


AJL 

0033 

FMPYP5 

AKL 

0035 

FN 

ANPJL 

0034 

FSUBD1 

ANPKL 

0036 

FSUBE1 

BINPNL 

0557 

FWDFFT 

BINPNT 

1176 

F5 

C 

0064 

H 

CLLRAR 

0667 

INDEX 

CPXMOV 

0741 

INDEXL 

CPXPML 

0554 

INXSTR 

CPXPMT 

0625 

J 

CTR1 

0760 

JLOOP 

CTR2 

0761 

JMSCPX 

Cl 3 

0415 

JMSFLT 

DATAST 

0020 

JMSFXP 

DATFLT 

0421 

K 

DJ 

1167 

L 

DNORM 

6600 

LFNOP 

DNORML 1164 

LLOOP 

DODFT 

0423 

LNOP 

D1 

0272 

LOOP1 

D5 

1364 

LOOP3 

ERRNN 

0471 

LOOP3L 

El 

0275 

LOOP4 

E5 

1367 

M 

FADDD1 

0445 

MAIN 

FADDE1 

0556 

MAINL 

FDIVFN 

0455 

MINJ4 

FENTER 

4407 

MINL 

FIXMOV 0701 

MINP 

FIXPML 

0416 

MOREL 

FIXPMT 

0600 

MP1 

FLIP 

0660 

MP1L 

FLPCTL 

0675 

M10 

FLPJ 

0677 

Ml 771 

FLPK 

0700 

N 

FLPLOP 

0665 

NMIN1 

FLTMOV 0721 

NN 

FLTPML 

10420 

PERMUT 

FLTPMT 

10606 

PI 


I 


0552 

PLOOP 

1006 

1361 

PMTDAT 

0614 

0555 

PMTMOV 

0653 

0446 

PMTNXT 

0644 

0402 

1 PRELIM 

0457 

1372 

Q 

1356 

0037 

1 REVFFT 

0524 

0262 

R4 

1173 

0024 

j R5 

1375 

0271 

1 s 

0067 

0031 

! SAVEJ 

0757 

1015 

: SETLOC 

0242 

0637 

’ SETLOL 

1355 

0613 

SETULL 

1354 

0605 

SETUPF 

1230 

0032 

SETUPL 

1073 

1171 

SPDFIX 

1166 

0553 

SPDFLT 

1165 

1013 

SPDNXT 

1136 

0624 

SPREAD 

1121 

1036 

SPREAL 

0417 

0200 

SWIT1 

0213 

1035 

SWIT1L 

0447 

1271 

SWIT2 

0222 

0025 

SWIT2L 

0450 

1000 

SWIT3 

1244 

0454 

SWIT3L 

0451 

1177 

SWIT4 

1250 

1172 

SWIT4L 

0452 

1170 

SWIT5 

1265 

1033 

SWIT5L 

0453 

0765 

T 

0072 

0560 

T3 

0762 

0522 

U 

0075 

0523 

UNLOOP 

1211 

0026 

UNSCRL 

0456 

0027 

UNSCRM 

1200 

0030 

ZP5 

1351 


064/ 

0021 



f 



28 






